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Abstract 

In this paper we consider the robust control of a ther- 
mal mixer using multivariable Sliding Mode Control 
(SMC). The mixer consists of a mixing chamber, hot 
and cold fluid valves, and an exit valve. The com- 
manded positions of the three valves are the available 
control inputs, while the controlled variables are total 
mass flow rate, chamber pressure and the density of 
the mixture inside the chamber. Unsteady thermody- 
namics and linear valve models are used in deriving a 
5th order nonlinear system with three inputs and three 
outputs. An SMC controller is designed to achieve ro- 
bust output tracking in the presence of unknown energy 
losses between the chamber and the environment. The 
usefulness of the technique is illustrated with a simula- 
tion. 


1 Introduction 

The problem of dosing the amounts of hot and cold 
fluid to produce a mixture having a desired temper- 
ature and flow rate is familiar to everyone, including 
people who are not control experts. In addition to the 
mundane shower example, the problem appears in in- 
dustrial practice as well [6, 9, 1]. In this paper we 
analyze the problem of obtaining the desired exit flow 
rate and temperature, while at the same time govern- 
ing the pressure of the mixture at a certain location. 
Of course, the inclusion of a third independent quantity 
to be controlled requires the addition of a third control 
variable, in this case a valve controlling the flow of mix- 
ture. Fig. 1 depicts the set-up. A mathematical model 
for the mixer in the context of propellant delivery in 
rocket engine test stands has been developed, and some 
controllers have been designed. In particular, small- 
signal linearization [1] and feedback linearization [9] 
controllers have been developed. In this paper we ex- 
pand the model to include first-order, linear dynamics 
for the valves and incorporate the effect of an unknown 
energy disturbance. In particular, we consider the dis- 
turbance to be the heat transfer between the chamber 
and the environment. This quantity is indeed difficult 
to measure and has practical relevance. Although pos- 
sible to model in principle, a quantitative description 
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of the heat transfer through mixer walls involves awk- 
ward convection correlations with parameters which 
are heavily dependent on temperature and flow charac- 
teristics. The convection coefficients can be extremely 
uncertain and often have to be found through iteration, 
even for the steady heat transfer calculations. For con- 
trol purposes, it is far more convenient to treat the 
heat transfer (and other unmodeled energy losses) as a 
lumped and unknown bounded disturbance. A bound 
on the external heat transfer is relatively easy to find. 
Other energy disturbances that intervene in the energy 
conservation equations can also be included to create 
a “bundled” uncertainty to be considered when tuning 
the controller. In order to simplify calculations, we as- 
sume that both cold and hot fluids are the same ideal 
gas obeying the well-known formula P = pRT , where 
P is the pressure, p is the density, T is the temper- 
ature and R is a constant particular to a given gas. 
As in [6, 9, 8], it is assumed that the valves conserve 
enthalpy and that source fluid properties, as well as 
exit pressure, remain constant . Although the temper- 
ature of the mixture past the exit valve is a controlled 
variable, it will be substituted by mixture density in- 
side the chamber to avoid cumbersome formulas. It 
will be assumed that the reference pressure and den- 
sity are chosen so that the desired exit temperature 
is obtained. This is possible since a given mixer pres- 
sure and density, along with exit pressure, uniquely de- 
termine an exit temperature through thermodynamic 
property data. Sliding Mode Control (SMC) is espe- 
cially suited for use in problems such as the one at 
hand, due to its insensitivity to certain kinds of distur- 
bances and parameter variations [3, 14, 13] and, in our 
case, the achievement of output decoupling, which facil- 
itates tuning. Applications of SMC to thermal systems 




have been reported, for instance, in [5, 11, 4, 12, 2]. 
2 Mathematical Model 


cold, hot and exit valves, and C v = [C vc C v h &ve ) T is 
the vector of valve flow coefficients. Using the above 
definitions and identities, the final form of the ideal 
gas mixer model with heat transfer disturbance can be 
written as 


As explained in [6, 9, 8], the transient versions of the 
mass and energy conservation equations are used in 
deriving the two basic model equations. It is conve- 
nient to write the equations in term of a pressure state 
instead of internal energy, since the flow rates can be 
directly computed from the state variables without fur- 
ther reference to thermodynamic relations. For the 
ideal gas, the following identities and definitions are 
in order 
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where h is the enthalpy, u is the internal energy, K p 
and K v are the constant pressure and constant volume 
heat capacities, respectively, and k is the heat capacity 
ratio K p /K v . As done in [7], we define a valve flow 
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F(P,p,C) = Cf(P,p) 


Thus, the nifctsa Sow nxies thiuiigh the two inlet valves 
and the exit valve are given by 


P = \c T v f (9) 

P = Ic^Apf+^Q ( 10 ) 

C v = A c C v + B c C (11) 

y = \P u C ve f e ] T (12) 


with A p given by Eq.( 13). 


A p = diag 
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and / = [/i / 2 — f e ) T - Q represents the rate of heat 
transfer (i.e., expressed in Watts in the SI system) be- 
tween mixer walls and the environment. 


3 Sliding Mode Control Design 

As seen in Eq. (10), the unknown heat transfer Q ap- 
pears in the equations as a disturbance to which the 
system can be made insensitive through SMC. Pres- 
sure losses due to friction can be combined with the 
heat transfer to create a bundled disturbance. For sim- 
plicity we choose three independent linear sliding man- 
ifolds, although other choices are possible [7]. Let the 
sliding manifolds be 


Fi (P, p, C vc ) = C vc fi ( P , p) for the cold fluid (5) 
F 2 (P,p,C vh ) = C vh f 2 (P,p) for the hot fluid (6) 
F e (P,p,C ve ) = C ve f e (P,p) for the exit fluid (7) 

where P and p are the pressure and density of the gas 
inside the chamber, and C vc , C v h and C ve are valve 
flow coefficients. Note that the input flows depend also 
on the thermodynamic properties of the source gases, 
regarded as constant parameters in Eqs. (5,6). The 
exit flow depends also on the pressure P s at the outlet 
of the exit valve, which is again considered as a con- 
stant parameter in Eq. (7). The flow formulas have 
a fairly general form, with the only requirement that 
they be linear in the valve position. We assume that 
the valves behave as a first order linear system between 
commanded position and flow coefficient. Valve dy- 
namics are expressed as 


si = Pd - P + c p (pd - p) (14) 

52 = P d - P + c P (P d - P) (15) 

53 = c w {w d — C ve f e ) (Ifi) 

where pd, P d and w d denote the desired trajectories 
for the density, pressure and exit mass flow outputs, 
respectively, and c p , cp and c w are sliding coefficients. 
Note that when the system is in the sliding regime, 
that is, Si = 0, i = 1, 2, 3, it is of degree five, the same 
as the uncontrolled plant. Therefore there is no zero 
dynamics with respect to the sliding manifolds and the 
system is trivially minimum-phase. Having established 
the stability of the zero dynamics, it is now possible 
to follow the steps of a standard MIMO SMC design. 
The first step is to write the 5th-order system in terms 
of output derivatives. Differentiating the first output 
twice leads to 


C v — A c C v + B c ( (8) 

where A c = diag^l/n, -l/r 2 , -l/r 3 ), B c = 

[/3 1 /? 2 /? 3 ] t , n are the time constants of the cold, 
hot and exit valves and (3i are their respective gains; 
( — [Ci C2 C3F is the vector of control inputs to the 
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Similarly, differentiating the second output twice gives 
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The third output, however, needs to be differentiated 
only once for the control to appear: 

2/3 
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It is desired that sj, s 2 and S3 reach zero after some 
finite time and remain zero despite of the heat transfer 
disturbance. Choosing Lyapunov functions 

K - is? 

for i = 1,2, 3, it is sufficient to ensure that V* < 0 at all 
times to guarantee the convergence of s*. To reach zero 
in finite time, a commonly used choice is to enforce 


Vi = SiSi < -T)i\Si 


(17) 


It can be shown [10] that the above implies that s t 
reaches zero in a time ti given by 


U < 


M* = Q)l 

Vi 


Enforcement of Eq. (17) can be achieved by using the 
equality for si and s 3 , since the heat transfer uncer- 
tainty is not involved in the expressions. Thus, let 
Si-ii = —r)i\si\ for i = 1,2. This is equivalent to 
Si = — r)isiga(si). Equality cannot be guaranteed for 
s 2 , since this would require knowledge of Q and Q. The 
SMC will rely, however, on the knowledge of a bound 
for these quantities and on the ability to increase the 
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ing Si = — r)isign(si) for i = 1, 3 leads to 
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where 
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and G is a 3-by-3 matrix with zeros in every entry, 
except for G(3, 3) = -c w (3z- For the remaining sliding 
variable, we choose 
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?2 = Pd + cpP d + r/2 sign (s 2 ) — 
f[((A. + =p/)3 p + ^)/ + A p f 


so that 


We want 
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with 7/2 > 0. Let 


BQ 

,a= ^ + m 

where Q is a known bound such that 

I Q + cpQ\ < Q 

at all times. Then it is straightforward to show that 
Eq. (21) is satisfied. Since 7/2 > 0, it is then sufficient 
to choose 

BQ , 00 , 

m > VK v ^ 


The equations that yield the control input, 
Eqs. (18), (20) and (19) can be compactly ex- 
pressed in matrix form as a system of equations linear 

in C 


(23) 


£_Bc 


' Ti ' 

f T A p B c 

V 

c = 

r 2 

f T G 


. r 3 . 


The existence of singularities preventing the inversion 
of the above system matrix and the more general ques- 
tion of state trajectory and control admissibility is cur- 
rently being studied by the authors in the broader 
context of constrained SMC. For Constant output set- 
points, however, it is possible to xinct conditions on 
the desired density, pressure and exit flow so that the 
steady state values of the valve flow coefficients are 
positive. The conditions reduce to [6] 

h c < h < hh 

where h is the enthalpy of the mixture that corresponds 
to the steady pressure and density setpoints, and h c 
and hh are the enthalpies of the cold and hot gases, re- 
spectively. Interestingly, this basic admissibility result 
is reduced to an expression involving only enthalpies. 
Note that the invertibility of the decoupling matrix for 
the feedback linearization design of [9] likewise reduces 
to 7— hh' 


Parameter 

Value 

Units 

Pi 

6.894 

MPa 

Ti 

311.1 

°K 

P 2 

5.516 

MPa 

t 2 

422.22 

°K 

Ps 

206.842 

kPa 

R 

296.8 

J/kg°Af 

K P 

1041.6 

J/kg °K 

k 

1.398 

n.a 

Po 

3.447 

MPa 

Po 

48.05 

kg/m 3 

CjU,2,3 

10 

n.a 

r l,2,3 

0.1 

sec 

01 

1 

n.a 

02 

2 

n.a 

fa 

3 

n.a 

V 

0.07079 

m 3 


Table 1: System Parameters for Sliding Mode Controller 
Simulation 


4 Simulation Example 

In this example, we choose the SI system of units and 
nitrogen as the working gas. The following design spec- 
ifications are used: 


• Obtain a steady exit flow of 10 kg/s with a mixer 
operating pressure of 4 MPa. 

• The temperature of the exit flow must be 
337 °K = 64°C'. 

• The heat transfer disturbance is Q = 10sin(10f). 

• The steady operating point must be reached in 
less than 25 seconds, for known initial conditions 
given in Table 1. 

• Valve chattering must not be observed. 


0 we obtain the reaching times: 

tr 1 
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The overall settling time is the sliding surface reach- 
ing time plus the settling time during sliding, given by 
t s i = t s 2 = t s 3 = 4/100 = 0.04sec.. The above gain 
choices thus satisfy the settling time requirements. Fi- 
nally, we check that the controller can attain its ob- 
jective despite the presence of the heat transfer distur- 
bance. The value of Q is computed as 

Q = \Q+c P Q\ = 1100 

We see that the selection of t ?2 is appropriate, since 
T& = 2.4062 x 10 6 > 

To avoid chattering, the signum functions present in 
the control law are replaced by finite slope approxima- 
tions consisting of saturation functions: 

sign(sj) » saX(si/<t>i) 


where 


(x if |aj| < 1 


sat(a:) = 


1 if x > 1 


1 otherwise 


(24) 


For the simulation example, the values <j>\ = <j>z = 10 
and d>2 = 1 x 10 4 were chosen. Figures 3, 4 and 2 show 
the simulation results. As it can be seen, the control 
objectives are satisfied. The second valve position is 
unacceptable, since it contains a time interval when the 
valve flow coefficient is negative. Avoidance of reverse 
flows and negative valve positions were not specified in 
the control problem. This simulation motivates further 
study into admissible state and control trajectories of 

— QA/rr> 

OVOtCillD UilU-Oi U1V1VJ. 


Other system parameters are summarized in Table 1. 

SMC Tuning The values of wj and Pd are directly 
taken from the design specifications. It is necessary 
to calculate the value of pd which will result in the 
desired exit temperature. For this, we have that the 
isoenthalpic process at the exit valve is also isothermal, 
since it is assumed that K p is constant. Therefore we 
choose pd = -jffr = 40kg/m 3 . In order to meet the 
settling time requirement, we choose c p = cp = c w = 
100, 771 = 773 = 120 and 7/2 = 2.4062 x 10 6 . With the 
known initial conditions and setting pd = pd = Pd = 


5 Conclusions 

The paper illustrates how MIMO SMC can be applied 
to the control of a thermal mixer in the presence of an 
unknown energy disturbance. The obtained controller 
guarantees robust tracking of the outputs as long as a 
decoupling matrix can be inverted. Invertibilty of the 
decoupling matrix is related to a broader problem of 
constrained control and is currently under investiga- 
tion. 
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Figure 2: Simulation of Sliding Mode Controller - Control Variables and Disturbance 

References 

[1] E. Barbieri and et. al. Small-signal point to point tracking of a propellant mixer. In American Control 
Conference , pages 2845-2850, 2003. 

[2] D. et.al. Das. Variable structure control strategy to automatic generation control of interconnected reheat 
thermal system. IEE Proceedings, part D: Control Theory and Applications , 138(6):579— 585, 1991. 

[3] U. Itkis. Control Systems of Variable Structure. John Wiley, 1977. 

[4] M.E. Jackson and Y. Shtessel. Sliding mode thermal control system for space station furnace facility. IEEE 
Transactions on Control System Technology, 6(5). 

[5] S. et. al. Ouenou-Gamo. A nonlinear controller of a turbocharged diesel engine using sliding mode. In IEEE 
International Conference on Control Applications, pages 803-805, 1997. 

[6] H. Richter. Modeling and control of E-l Hydrogen mixer; Model development. Technical report, NRC Research 
Associate, NASA Stennis Space Center, Mississippi, 2002. 

[7] H. Richter. Tracking of a thermodynamic process using a polytropic surface as sliding manifold. In American 
Control Conference, pages 197-201, 2003. 

[8] H. Richter and et. al. Modelling, simulation and control of a propellant mixer. In AIAA/ASME/SAE/ASEE 
Joint Propulsion Conference and Exhibit, 2003. 

[9] H. Richter and et. al. Nonlinear modeling and control of a propellant mixer. In American Control Conference, 
pages 2839-2844, 2003. 

[10] J.J. Slotine and W. Li. Applied Nonlinear Control. Prentice-Hall, 1991. 

[11] M. et. al. Sommerville. On variable structure control of a throttle actuator for speed control applications. In 
IEEE Workshop on Variable Structure Systems, pages 187-192, 1996. 

[12] S. Spurgeon and Ch. Edwards. Nonlinear control of a 175 kw gas burner. Computing and Control Engineering 
Journal, 6(l):29-34, 1995. 

[13] V.I. Utkin. Variable structure systems with sliding modes: A survey. IEEE Transactions on Automatic 
Control, AC-22:212-222, 1977. 

[14] V.I. Utkin. Sliding Modes in Control Optimization. Springer- Verlag, 1992. 







Sliding Functions 



Time, sec. 

Figure 3: Simulation of Sliding Mode Controller - Sliding Variables 
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Figure 4: Simulation of Sliding Mode Controller - Output Variables 
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